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Summary 


The following is a list of our accomplishments in the past 12 
months. Most of the objectives in our original proposal for this period 
have been achieved. Additional research was also performed at the 
request of NASA/GSFC, code 924. 

1. Receiver performance study of space-borne laser altimeters and 
cloud and aerosol lidars: 

We have developed a procedure about how to simulate the APD 
output signal using the much more exact models given by McIntyre and 
Webb for the APD. The preamplifier noise model was also improved. This 
APD simulator will be used to update the existing Laser Altimeter 
Simulator which currently uses an over simplified Gaussian model for 
the APD output photocurrent. A separate write-up titled "Computer 
Simulation of the output signal from an APD" is attached. 

A computer program was written to simulate the noise counter 
output of the laser altimeter. The nearly exact models for the APD and 
preamplifier noise were used. The noise counter output serves as the 
input to the automatic threshold adjustment circuit of the laser 
altimeter. This program is currently being incorporated in the Laser 
Altimeter Simulator, version 4.0. A short report about our program 
titled "Simulation of the number of threshold crossings due to noise 
between laser shots for the GLAS simulator" is included in this report. 

A computer program was written using the software Mathcad to 
analytically calculate the performance of APD laser altimeter receivers 
based on Webb's distribution for the APD output. The Gaussian model 
was shown to have series shortcomings. The program used the 
Geoscience Laser Altimeter System (GLAS) as an example. A copy of the 
program, "GLAS altimeter link margin analysis", is attached. 

A simple laser ranging system was set up in the lab. A low power 
1.06 (im wavelength laser diode was used as the test laser transmitter. 

The use of a laser diode allows us to arbitrarily alter the laser pulse 
shape in order to simulate the altimeter return signals. The total timing 
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jitter of the pulse generator, the laser, and the counter/timer under 
strong input signal conditions was found to be about 60 ps, which was 
well within 300 ps required by GLAS. 

We conducted with D. Reusser of USRA a study of radiation damage 
to Si APDs for the radiation environment similar to that of GLAS. We 
found that all previous radiation damage tests used either electron or 
gamma rays as the radiation sources. The radiation of concern for low 
earth orbit such as GLAS is primarily high energy protons at the so 
called South Atlantic Anomaly region and poles of the 
geomagnetosphere. Based on our study, Si APDs were expected to suffer 
long term radiation damage, which may be negligible for linear mode 
operation but overwhelming for Geiger mode operation. A proton 
radiation damage test for silicon APDs was planned. 

We also performed a research of available Si APDs as candidates 
for the GLAS altimeter receiver. A custom APD preamplifier module by 
EG&G was chosen. Three devices were purchased and they were 
scheduled for delivery in May 1994. 

The receiver signal to noise ratio for the GLAS cloud lidar was 
calculated assuming a PMT photon counter is used as the photodetector. 
The total charges drawn from the PMT photocathode and anode were also 
computed. A copy of the Mathcad program printout, "GLAS cloud lidar 
receiver photon rate calculation", is included in this progress report. 

2. Receiver performance analysis for space-to-space laser ranging 
systems: 

All our work in this area was related to the Gravity And Magnetic 
Experiment Survey (GAMES) mission proposed by NASA/GSFC, code 924. 

We studied different types of receiver structures for GAMES. The 
most attractive receiver design appeared to be the one used by the 
National Radio Astronomy Organization (NRAO), which used a down- 
converter and a narrow band filter instead of a PLL. NRAO had 
successfully used this type of receiver for laser ranging to stationary 
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objects. The receiver can be greatly simplified when using this approach 
as compared to other types of receivers. 

We studied the receiver performance for a moving target. The 
receiver estimation bias due to Doppler shift was analyzed. The effects 
of the finite sample rate and number of bits of the analog to digital 
converter were calculated. The results were consistent with the actual 
measurements. A short report on this subject, "Use of the NRAO method 
for the GAMES fine ranging receiver 1 ', is attached. 

The link analysis of the GAMES fine ranging system was improved 
and included the use of APDs. The effects of the PMT gain and stray light 
were calculated. The receiver performance using a PMT and an APD was 
compared. The program was used to calculate the performance of the 
latest version of the GAMES laser ranging instrument (LRI-2). 

The receiver design for a multitone laser ranging system was 
studied. A copy of our results, "Receiver performance analysis of a 
multiple tone laser ranging system", is attached. 

CCD cameras, star trackers, and quad APDs as the tracking 
detectors for GAMES were also studied. 

We have estimated (with D. Reusser, N. Walsh, and K. Mehalick) the 
amount of stray light incident on the fine ranging and coarse ranging 
detectors. The major stray light source was found to be sunlit clouds. 

The sun glint from the ocean was shown to be negligible. The moon and 
other planets and stars were found to have little contribution to the 
total stray light onto the detector unless they are in the telescope field 
of view. The minimum angle was found between the sun and the 
telescope line of sight for the receiver to remain operational. The 
details of this study were already reported in our progress report for 
the previous six month period. 

A preliminary study of radiation damage to the PMT and CCDs for 
GAMES was completed (with D. Reusser). The major radiation damage to 
PMTs occurred on the glass window. PMTs have been tested and shown to 
suffer little performance degradation for the radiation dose experienced 
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by GAMES. CCDs were found to be relatively sensitive to radiation 
damage. Test data in the literature showed that regular CCD cameras 
would suffer serious long term damage if used in GAMES and the 
resultant performance degradation would not be acceptable. A well 
designed rad-hard CCD camera, on the other hand, has been shown to 
suffer negligible long term damage for the radiation dose experienced 
by GAMES. A separate report, "Space radiation effects on photodetectors 
for GAMES and GLAS", is attached, which included our results on PMTs, 
CCDs and Si APDs. 

Pseudo noise (PN) code laser ranging systems were studied for 
the GAMES coarse ranging subsystem. Receiver performance was 
analyzed for both baseband and subcarrier modulated PN ranging codes. 
The so called JPL ranging codes were studied as a possible candidate for 
the GAMES coarse ranging system. The power spectra of regular PN 
ranging codes and JPL ranging codes were calculated. Included in this 
report are copies of two short reports, "Receiver performance analysis 
of the GAMES PN code coarse ranging subsystems" and "Power spectrum 
of PN ranging codes." 

3. Receiver performance study for the Mars Environmental Survey 
(MESUR): 

We have done preliminary estimates of the solar background 
photon flux. It is expected that most of the detected photons will be of 
solar origin. The number of bits required to store each histogram bit and 
the signal to noise ratio both depend strongly on the number of detected 
signal photons, the solar photon flux, and the averaging time. Therefore 
specification of the histogram electronics and averaging time depend on 
an accurate estimate of the solar flux. 

We have investigated several low power, radiation tolerant 
computers. A Zilog Z80 has be selected for prototyping. In addition, a 
low power memory has been selected. 

Different types of quenching circuitry for the Geiger mode APD 
photon counters were studied and tested for the lowest electrical 
power consumption. 
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Computer Simulation of 
the Signal Output from an APD 

Xiaoli Sun 

The Johns Hopkins University 


1. Introduction 

A computer simulator of space borne laser altimeters which use 
silicon APD has been developed by NASA/GSFC. However, the signal 
output from the APD was simulated approximately using the 
Gaussian distribution. The actual distribution of the signal output 
from an APD was given by McIntyre [1] which was significantly 
different from the Gaussian distribution, especially at the tails of the 
probability density functions. A close approximation of the McIntyre 
distribution has been given by Webb [2]. This report attempts to give 
a guide line of how to upgrade the existing laser altimeter simulator 
with the much more exact APD models by McIntyre or Webb. 

The signal to be simulated consists of photoelectrons output from the 
APD and noise from the preamplifier. The former follows the 
McIntyre distribution and the latter follows the Gaussian 
distribution. The two are independent of each other. The output of 
the simulation program consists of discrete samples at a specified 
sampling rate. Each sample is equal to the signal averaged over the 
sampling interval. From now on, we define the signal as the input 
current to the preamplifier. All the noise sources will also be 
converted to their equivalent noise currents at the input to the 
preamplifier. 
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2. Simulation of the Preamplifier Noise 

The preamplifier noise follows the Gaussian distribution. Press et al 
[3] have shown subroutine for generating standard Gaussian random 
variable. We need only to scale this standard Gaussian random 
variable, y g o, to obtain the preamplifier noise output, y g , according to 
the relationship 

y t = °y t0 +Li (i) 


where p and a are the mean and standard deviation. The mean of 
amplifier noise can be considered to be zero. The variance of 
preamplifier noise can be calculated as 

, 1 r- sin 2 (Q)Ar/2) ... 

cr = — T-^Ni^da) 

27T j 0 (coAt/2) 2 


( 2 ) 


where N(co) is the one sided noise power spectrum of the 
preamplifier and Ax is the sampling interval. The first part of the 
integrand in (2) results from the fact that each sample is the average 
signal over the sampling interval (integration of the actual signal 
over Ax and then divided by Ax). If the preamplifier contains a FET 

as the front stage, the noise power spectrum can be written as [4] 
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( 10 ) 


where K is the Botzmann's constant, T a (°K) is the ambient 
temperature, Rf is the transimpedance of the preamplifier, g m is the 
transconductance of the FET, T is a factor close to unity, and Ci is the 
total capacitance at the input of the preamplifier. 
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When the altimeter receiver bandwidth is below 200 MHz, one can 
often approximate the preamplifier noise as white Gaussian noise 
with the variance given as 


0 * = 


4AT. 1 

Rj 2Ar 


( 11 ) 


where T n (°K) is the equivalent noise temperature. 

3. Simulation of the Photocurrent Output from the APD 


The sampled output of an APD can be modeled as a discrete random 
process which is equal to the number of photoelectrons during the 
sampling interval times the electron charge and then divided by the 
sampling interval. The probability of the number of photoelectrons 
has been shown to follow the so-called McIntyre distribution [1], as 
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t?i n\ 

/ > (0IA(/)At) = e" a(,)AT (12) 


with Pm (min) the probability of generating m secondary 
photoelectrons in response to n primary photoelectrons, as 
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where k c ff is the APD hole to electron ionization coefficient ratio, G is 
the average APD gain, and X(t) is the primary photon counting rate 
given by 

Ut) = £[P,(t) + Po] + ^ (14) 

hj q 


where rj is the APD quantum efficiency, hf is the photon energy, P s (t) 
is the input optical signal power, Po is the background radiation 
power assuming to be constant, lb is the APD bulk leakage current, 
and q is the electron charge. 


A very close approximation to (12) was given by Webb et al [2] for 
A(t)Ax»l, as 
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(15) 


where F is the APD excess noise factor which can be calculated as 


F = ^G + 



(16) 


It also known [5] that the Webb’s approximation becomes inaccurate 
as A(t)Ax»l. Therefore, it is suggested to use Eq. (15) only for 
A(t)Ax >10. 


Because the complexity of the probability distributions for the APD, 
one cannot generate the random variable as an analytical function of 
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a uniform distributed random variable. Instead, the rejection method 
described in [3] has to be used. A comparison function has to be 
found when using the rejection method. The comparison function is a 
function which is greater than the desired probability distribution 
function for all possible value of the random variable and the total 
area under the curve is finite. The inverse function of the indefinite 
integral of the comparison function should also be known 
analytically. The area (integral) of the comparison function should 
also be as small as possible since the time required to generate the 
random number increases with the area of the comparison function. 

We may simulate the APD photoelectron output for X(t)Ax < 10 in 
two steps. First we generate a Poisson variables to simulate the 
number of primary photo electrons. We then simulate how many 
secondary photoelectrons each primary photoelectron will generate 
according to Eq. (13) with n=l. The APD output should be the sum of 
secondary photoelectrons for all primary photoelectrons. This 
procedure involves calling a Poisson random number generator once 
and then calling a McIntyre random number generator n times with 
n the outcome of the Poisson random number generator. This method 
only requires a comparison function for the McIntyre distribution for 
the n=l case. 

The Poisson random number generator can be implemented using 
the computer subroutine provided by Press et al [3]. The McIntyre 
random number generator can be realized similarly but with the 
propitiate probability distribution and comparison functions. We may 
still use a Lorentzian distribution as the comparison function, i.e., 
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/(*) = 
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(17) 


We need to find a set of values for co, xo, and ao such that 
f(m)>PM(mll) and J/(x)dr as small as possible for all possible values 

of m. One example is co=4/G, xo=-G/50, and ao=G/3. Appendix A 
shows plots of Pm (mil) and f(x) with the above parameter values for 
various APD gains from G=50 to 500, 

For X(t)At >10, one may simulate the APD output in one step using 
the Webb's approximation (15) to save computer times, A suggestion 
for the comparison function is still the same as Eq. (17) but with 
co=l/2G(FX(t)AT)l/2 # xo=GX(t)At, and ao=G(2FX(t)At)l/2 > Appendix B 
shows plots of the Webb's approximation and the suggested 
comparison function for X(t)Ai =100 and various values of the 
average APD gain. 


APDs also have surface leakage current which is not multiplied by 
the APD gain. This part of the APD dark current should be modeled 
as a DC current with its shot noise contribution following the 
Gaussian distribution with the mean and variance given by 


H = ql. At 
of = qI,Ax 


(18) 


where I s is the surface leakage current. This noise can be added to 
the total preamplifier noise. 
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Simulation of the Number of Threshold Crossings 
due to Noise between Laser Shots 
for the GLAS Simulator 


Xiaoli Sun 
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January 1994 


1. Theory 

The threshold crossings due to noise, or false alarms, can be assumed 
to follow the Poisson distribution for a reasonably low false alarm 
probability. The rate of false alarms can be calculated as [1] 

T - 1 -P fa ( i \ 

1 FA T~ ( 1 ) 

where pf a is the probability of threshold crossing at a fixed time 
within the observation interval due to noise alone and x is the 
average pulsewidth of the false alarms, which is roughly equal to the 
reciprocal of the receiver bandwidth. The total number of false 
alarms in an observation time interval, T, can be written as 

P P A(n)=^fLexp{-N FA ). (2) 

where Nfa is the average number of false alarms given as 

^FA = ^^ = TPfa (3) 

The point false alarm probability, pf a , can be calculated as 

p fa =[ p(xlO)dx (4) 

JXj 

with p(xlO) the probability density function (p.d.f.) when no signal is 
transmitted. For convenience, we define the random variable x as the 
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equivalent number of photoelectrons integrated over x second in 
front of the preamplifier. 

\ 

There are two independent noise sources, the APD dark noise and the 
preamplifier noise. The p.d.f. of the total noise is given as the 
convolution of the p.d.f. of the two noise sources. 


The number of equivalent dark photoelectrons from the APD has 
been shown to follow the Mclntyre-Conradi distribution [2] [3], which 
is rather complicated. A much simpler and yet very close 
approximation to the Mclntvre-Conradi distribution has been oiven 
by Webb et al. [4], as 


P pd( x !Po) = 
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expJ 
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(5) 


where G is the average APD gain, F is the APD excess noise factor, 
and po is the average number of equivalent primary APD dark 
photoelectrons given as 



( 6 ) 


where nb is the average background radiation noise count over a 
laser pulse width time, I b is the APD bulk leakage current, and q is 
the electron charge. The random variable, x, in (5) should be a 
discrete number. But the average APD gain is often set to a relatively 
large value for high receiver sensitivity and the distribution of x 
may considered as continuous. The APD excess noise factor in (5) can 
be calculated as [4] 


2 


(7) 


F=G.k eff+ (2-l)(l-k ef( ) 

where k c ff is the APD ionization coefficient ratio of holes and 
electrons. 


The noise contribution from the preamplifier appears additive and 
follows the Gaussian distribution with zero mean and variance given 
by [5] 


2KT n x I, 

R,q 2 q 


where K is Boatsman constant, T n (°K) is the equivalent noise 
temperature of the preamplifier, and Rj is the APD load resistance 
(the preamplifier gain). We have assumed here that the amplifier 
noise is white noise, which is a good approximation up to a 
bandwidth of about 200 MHz. The APD surface leakage current also 
contribute additive white Gaussian noise and it can be included in 
the equivalent noise temperature. 


The p.d.f. the total noise is the convolution integral of Eq. (5) and the 
p.d.f. of the Gaussian noise, i.e., 
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The probability of false alarm can be obtained by substituting (5) 
into (9) and then (4), as 
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The double integrals in (10) can be simplified by first interchanging 
the orders of the integrals and then using the built-in error function 
of the computer. That is. 


Pfa= 


/- W X T 


J ^2<j^ ^ X | P pD (y^o)dy-J Ppo(yipo) *(V 


where 




Using the built-in complimentary error function, 
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Substituting (13) into (11) 


Pfa=| PpoCyipo) ^ erfc (^^jdy 


X )dy (11) 
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2. The program 


Figure 1 shows a block diagram of the Fortran program which first 
evaluate Eq. (14) and then generate a Poisson random number with 
the mean given by Eq. (1). The subroutines and functions used for 
the numerical integration and Poisson random number generation 
are copied from those given in [6, ch. 4 and ch. 7], 
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Figure 1 . Block diagram of the program. 
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Calculation of the detected signal and noise photons per pulse 
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Other system parameter values: 
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The p.d.f. of the APD output: 
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NT,, 


5 

5.5 

6 

6.5 

7 

7.5 

8 

8.5 

9 

9.5 

10 


P^k 

3.3545- 10" 5 
8.6669- 10~ 6 
2.1431- 10~ 6 
5.0946- IO" 7 
1.1689- 10~ 7 
2.5971- 10~ 8 
5.6053- 10" 9 
1.1783- 10~ 9 
2.4181- lp" 10 
4,8545- 10" 11 
9.5513-10" 12 


PFA k 

0.2004 

0.0561 

0.0142 

3.3906- 10~ 3 
7,7894- 10~ 4 
1.7313- 10~ 4 
3.7368- 10~ 5 
7.8552- 10~ 6 
! 1.612- 10~ 6 
3.2364- 10~ 7 
6.3675-10“ 8 


TOL •'= 10~ 7 

Rj N^Rj) Pp(N s ( Rj), 7.0) 


2.5- 10 2 

3.7393- 10 3 



tHEM 






5-10 2 

9.3483- 10 2 

5.5- 10 2 

7.7259- 10 2 

610 2 

6.4919- 10 2 

6.5- 10 2 

5.5315- 10 2 

7-10 2 

4.7695- 10 2 

7.5- 10 2 

4.1548- 10 2 

8-10 2 

3.6517- 10 2 

8.5- 10 2 

3.2347- 10 2 

9-10 2 

2.8853- 10 2 

9.5- 10 2 

2.5896- 10 2 

1-10 3 

2.3371- 10 2 

1.05- 10 3 

2.1198- 10 2 

1.1- 10 3 

1.9315- 10 2 

1.15- 10 3 

1.7672- 10 2 

1.2- 10 3 

1.623- 10 2 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

0.9997 

0.9964 

0.9771 

0.9183 

0.8058 

0.6524 

0.4893 













GLAS Cloud Lidar Receiver Photon Count Rate Calculation 

Xiaoli Sun, Johns Hopkins University, March 1994 


Instrument parameters: 

Laser pulse energy (Joule/pulse): 

Laser wavelength (nm): 

Laser pulsewidth (s): 

Laser pulse repetition rate (Hz): 

Laser beam divergence angle (rad): 

Lidar altitude (m): 

Receiver optics FOV (rad): 

Receiver telescope diameter (m): 

Telescope secondary mirror diameter (m): 

Total transmission of the receiver optics: 

Receiver optical bandwidth (nm): 

Range gate (km): 

PMT quantum efficiency: 

PMT gain: 

PMT dark counts/s: 

Receiver photon counting interval (s): 


electron charge (C): 
Speed of light (m/s): 

Photon energy (Joule): 
Range gate interval (s): 


X 


RG- 10 

T rG 2 


E t := 0.050 
X := 532 

x :=20 io " 9 

Fp : =40 

e d :=iooio ” 6 

R := 705- 10 3 
0 pov " = 250- 10 ^ 
♦ tel := °- 90 
4>obs := 0- 2 0 
T 0 :=0.5 
AX :=0.10 
RG := 20 

t\ := 0.20 
g := io 5 

nd pmt * = 150 
At := 0.5-10" 6 

ql= 1.6-10” 19 
c := 3 - io 8 

T RG = 1.333- 10 -4 



Atmosphere parameter values at 532 nm wavelength: 


Solar spectral irradiance (W/m A 2 um): 
Lunar Spectral irradiance (W/m A 2 um): 

Albedo of the terrain: 

Albedo of the cloud: 


E X solar := 1848 
E Xlunar := 00044 


1 ter :=0 - 34 


r cloud 


:=0.70 


Spectral radiance intensity of terrian at daytime (uW/cm A 2 um sr): 


day 


l-'-r 


'dhy' 2-10 


Spectral radiance intensity of sunlit cloud (uW/cm A 2 um sr): 

io 2 ) .LSpS 


I a — i — a \sQjaf - 1( 


I daycloud 4- ^ 8 


daycloud 

Maximum spectral radiance intensity at night at full moon (uW/cm A 2 um sr): 


1 night = °-°^ 8 


Aerosol cross section times the square of the atmosphere transmission (1/km sr) 
as a function of altitude (km): 

i := 0, 1 .. 48 HT := READPRN( GLASHT) pT2 := READPRN( GLASBTT) 

pTT(h) :=linterp(HT.pT2,h) h :=0.1,0.2.. 25 
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Calculated Results 


Received signal photons/s: 


Received laser signal power (Watts): 
71 • (<t> let 2 ” <t> obs 2 ) 


p «(h) : =T E f 


4-(R- h-lOOOr 


T 0 (prr(h)io" 3 ) 


Peak signal (Watts) (at 9.5 km altitude): 


P j(9.5) = 5.1572* 10 


“11 


Detecteded signal photons/s as a function of the cloud altitude in km: 


n .(h) :=tv 


P r (h) 


hf 


n _( 9.5 ) = 2.761*10 


Detected signal photons/s as a function of time: 

Reflection from the ground: 

E t r ter n ' ($ tel 2 “ obs 2 ) 


nl Q (t) :=if1 


/ R \ . b t r ter 

t> 2— + x ,0, 

\ c / x n 


4R^ 


Total signal photons/s: 


0.5T 0 .o(,-2.|).i 
nl s ^2 — + 10" j = 4.407*10 


ns( t) :=if|t>^2^ + 2010' 9 j,0,n s ^R--yj l0" 3 +nl s (t)J 

Background radiation noise photons/s: 


Receiver FOV solid angle: 

Day time: 

Clear Sky: 


n := 


«• 0 FOV 


Pnd-Oday 10 " 6 ) ■■■ — * 0bS ^ 10 T 0 n (AXlO~ 3 ) 


K - PBd 

nb d .-n — 


PB d = 2.969* 10 


nb d = 1.589*10 


-10 


Bright cloud: 


-6\ " (‘frtel 2 ♦ obs 2 ) ' 1q4 orAJLin" 3 ! 
dcloud ) 4 Tq Q (AX 10 ) 


PB 10 


nb 


dcloud ■"* T 1‘ 


P® dcloud 


PB dcloud = 61 12 * 10 


“10 


hf 


nb dcloud = 3272 * 10 




Maximum PMT current (A): 

Cathod: I pMTc := ( nb ±kMd + n s ( 9.5 ) )• q I pMTc = 5.678- 10' 

Anode: 1 PMTa := 1 PMTc G 1 PMTa = 5 678 ‘ 10 


Upper bonds for the total electron charge drawn from the 
PMT photocathod: 

Average number of detected photons 
per 24 hours with the range gate: 


n := 

n s ( 9.5 ) + 0.5- 

i n i 

+ 0.5- nb n 



i 2 j 



•TRG-Fp-3600- 24 


N = 3.604* 10 


10 


If there were always bright clouds: 


N c := 

n s ( 9.5 ) + 0.5- 

dcloud' 

+ 0.5- nb n 



1 2 i 



TRG-Fp-3600- 24 


Average number of detected photons per day with NO range gate: 

•3600-24 


N c = 6.072* 10 


10 




f nb dl 


n nrg := 

n s ( 9.5 ) + 0.5- 

ill 

+ 0.5- nb n 


n nrg = 6 - 757 * 10 


12 


If there were always bright clouds: 


N NRGcloud := 


n.( 9.5 ) + 0.5- 


/ nb dcloud' 

+ 0.5- nb„ 

« 

\ 2 J 



3600-24 


N NRGcloud " U39 ' 10 


Total number detected photons over the GLAS 5 year mission: 

With range gate: 

N 5yr :=N c -5-365 N 5yr = 1.108* 10 14 

With NO range gate: 

N NRGSyr := N NRGcloud' 5 ' 365 N NRG5yr = 2 - 078 ‘ 1q1 6 

Total electron charges drawn from the PMT photocathod (Coulomb): 


With range gate: 

C cathod • =N 5yr' 1 - 6 ' 10 19 c = 1.773* 10 _5 

Without range gate: 

C NRGcaihod := N NRG5yr‘ 1 -6- 10 1 9 C NRGcathod = 0 00332 
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Use of the NRAO Methodt 
for the GAMES Fine Ranging Receiver 

Xiaoli Sun 

The Johns Hopkins University 
Baltimore, Maryland 21218-2686 

February 1994 

1. Analysis of Phase Estimation Bias due to Donpler Shift . 

The average received analog signal from an ideal photodetector may be 
written as 


*.(/) = 4 cos[(a) e + co D )t + 6] 


( 1 ) 


The actual signal into the computer is the truncated and digitized signal 
given in ( 1) and it can be written as 


N+riQ 


f 


\ 


x(t)~ x n 8(t-nM) = w(t )■ ^S(t-nAt) 


xjt) 




\n=— • 


( 2 ) 


where 


w(t) = 



t 0 <t<T + t 0 
Otherwise 



(3) 


and 


x n = x 0 (nAt) = Acos[(co c + co D )nAt + 6] 


( 4 ) 


t J. M. Payne, D. Parker, and R. F. Bradley (National Radio Astronomy Observatory, 
Charlottesville, VA 22903), "Range finder with fast multiple range capability," 
Rev. Sci. Instrum., Vol. 63, No. 6, pp. 331 1-3316, June 1992. 
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The phase esitmator using the NRAO method is given by 


g =tan -[MM, 


> = tan 


-i 


sin(o) c nAr) 

« = n 0 

A+* # 

^x k cos( co c n At) 


(5) 


where 


X(co c ) = ]x{t)e' ia ‘'dt 

= ^ x n Je” ; “’‘'<5(r - nAt)dt (6) 

— 

N+n 0 A + a, 

= ^x H cos(co c nAt) + ; sm(£y c nA/) 

« = r =« # 


Eq. (6) should also be equal to the Fourier transform of the right hand 
side of Eq. (2) which can be written as the convolution of the Fourier 
transforms of the rectangular window function, the sample function, and 
the signal, as 


X(co) = W(co)* 




k=-° 


(7) 


Since 


and 


W( co) = T ■ sine 


. fcoT' 


{■ 


y-J^o 


X 0 (co) = -{e J9 8[co-(co c +co D )] + e~ j9 6[co + (co c +co D )]}, 


(8) 

(9) 


therefore 
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X(o» = 


kTA 

At 


, , 2tz 

- (CO-k T .2a. 

X , At , o 

Sl nc( — Je 


Al 




*{e je S[ co-(o) c + co D )J + e- je 8[ co + (co c + co D )]} (10) 


X(co) = 


kTA 

At 


e jd ^ sine 

* 


2 TC 

(co-(co c +co D )-k—)T 
At 


-j(u-(a c + a D )-k.llL)t a 
e ** 


+ e iB ^ sin c 

k=-°° 

2k 

(CO + (CQ c + CD D )-k )T 

At 

> 

-J( a-{co c -r )-k^A)t 0 

e to 

2 


- 

> 


( 11 ) 


2k 


Assuming — - » co c + co D and filtering out all the k* 0 terms, 


At 

X(co)- 


^\e Je sinc(ttsL) e '“«* 

Ar L ~ 


+ £ 


(26) c + (DqJT- 
2 


,-;Y2<a e +<u i) ;;/ 0 1 


( 12 ) 


The phase esitmator becomes 


e = tan l * 

Im[X(co e )f 



[Re(X(a> c )J J 


= tan~ ] < 

r 

sin((O D t 0 + 6)sinc( D ) sin[(2co c + co D )t 0 + 9] sine 

\(2<o e +<D B yr~l 
2 


co T 

cos(co D t 0 +9)sinc( D ) + eos[(2o) c + co D )t 0 + G] sine 

l j 

(2co c +co D ) rl 
2 


( 13 ) 

If one select t 0 such that 2 co c t 0 = 2mK with any integer m, then 
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sin[9 + (2co c + co D )t 0 J = sin(9 + co D t 0 ) 


(14) 


and 


9 = tan 


sin( co D t 0 + 9) sin c (— ) - sin( 0) D t 0 + 9) sin c 

~(2(D c + cq d )T~ 
2 

cos(co D t 0 + 9 )sinc(—^~) + cos( co D t 0 + 9) sine 

~(2cd c + co d )T~ 

2 J 


= tan 1 < 


, co d T / 2 

tan(co D t 0 + 9) - 


sin(co D T/2) sin[(2co c + co D )T / 2] ] 

(70),. + Ct)rJT /2 


j/nf fit) D r / 2,) j/n/~ (2 <y c + JT 7 27 

0^1/2 (2co c + cd d )T 12 J 


Tjr ^ „ . . r (2co r + co D )T 1 . ,co n T . 

If (D c T -2nK for any integer n, sin[ — ] = sin( — — J and 

2 2 


(15) 


B — tan 


-1 


1 


1 


= tan"'<( 


tan(co D t 0 + 9 ) 


tan(<v o + 0) 


(2a) c -+&),,) 

1 1 

1 , 

© X , (2co c + co d ) J 


fit). 


a,+fit) 0 


(16) 


or 


tan 0 = tan(fit) c r 0 + 0)' 


<y. 


0). + £t). 


(17) 


A 

The relative error in tan 9 is given by 


e- = 


tan# -tan (fit),/. + 0) co 


^0 _ ^ D 


tan (o^ o +0) 


<y + fir 

c 


fit) 


(fit)a«0) c ) (18) 


Since d tan 9 = r , the phase estimation error can be written as 

cos 9 
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Afl - cos 0 . A tan § = cos 9 • tan(*y 0 4- 0 ) . + 

tan(fi) D / 0 + 0) 

= cos(ov 0 + 0) • tan(GVo + 0) — — — 

^c+^Z, 

co n 

^-^■sm(co D t 0 + 6) 

The relative phase error is given by 


AO 


co D sm{co D t 0 +9) 



co D t 0 + e 


v D t 0 + o 

X 

J 

CO 

C 


(19) 


( 20 ) 


An unbiased phase estimator is suggested according to (16), as 

X*. sinCo^nA/) 


0 0 = tan" 1 


Uc±Uj> 


CO 


A’ 4- 

X-r, cos(o) e nAr) 


(21) 


where is the estimated Doppler shift obtained from the previous 
measurement. 

2, Minimum Required Sampling Rare 


The sampling rate after the down converter has to be fast enough to 
avoid aliasing. The Fourier transform of the sampled signal is given by 

Eq. (11) and it is plotted in Figure 1. It is obvious that the sampling rate 
should be such that 

2 x 0 / s 2k 

— -2(ao e + (Q D )»— (22) 

For GAMES, the carrier frequency is 10 KHz (co e = 2;rxl0 4 ) and, the 
maximum Doppler shift is a few Hertz, and the observation time is T=0.1s 
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a sample rare of 3 times the carrier frequency, or 30 Ks/s, will be 
sufficient. 


3. Quantization Caused bv the A-to-D Converte r 


Quantization errors result from the limited resolution of the A-to-D 
converter. In this section, we analyze the maximum quantization error 
assuming there is no Doppler shift. 


Considering the quantization error, the phase estimator given in Eq. (5) 
should be rewritten as 




6' = tan' 1 


Z(*. + Ax K )sin((O c nAt) 

A+n # 

+ AxJcos(co e nAt) 




A’ 


= tan 




^x n sin(co c nAt) + Ax, sin(6J £ /iAt) 

A T *f« # A + * # 

^x H cos(co e nAt)+ ^Ax K cos(co c nAt) 


(23) 


where Ax n 's are the quantization errors. One needs only to consider the 

^ A* + A t +h # 

case when 6' = 0. Under this condition,^ x K sm(co c nAt) « ^x n cos(co e nAt) 


n= 




A r +% AV« # 

^x K sin(co e nAt) ^Ax K sm(co c nAt) 

/ *=% r.-n* 

cos(fii f nA/) cos(a) e nA/) 


(24) 


If the resolution of the A-to-D converter is given as Ax, the maximum 
uncertainty of each sample is 0.5Ax, or, lAt n i £ Ax/2, the phase error due 
to quantization error becomes 
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A 0 = \9'-e 


N*r, 

sin(co c nAt) 


N+* m 


2 x „ cos(o) c nAr) 

A*. S| sin ( £y c /2A/ )| 


A'+n 


^ r B 


A’+n^ 

X*, C0S(&) c /7A/)j 


^£.w. 2.-1 

2 &j e __ 


.4 


r*>. 


J cos(o) c / + 9)cos{co e t)d : 


2Ax 


A- — -cos(9)-T xAcvs(d) 
2 


_2 Ax 
k A 


(25) 

where the observation interval, T, is set to T=N>(27i/co c ). For a given 
maximum allowed phase error, A 9 mtx (rad), the resolution of the A-to-D 
converter has to satisfy 

Ax<jA-A9 wil (26) 


The minimum number of bits required for the A-to-D converter is 


N, 


bits 



> lo<* 



(27) 


For example, a 9 bits A-to-D converter will give a maximum phase 
quantization error of 3.1 mrad, or 1/2048 of a cycle, which is sufficient 
for GAMES (50 pm at 2 GHz which corresponds to 4.2 mrad orl/1500 
cycle). 
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November 10, 1993 


X. Sun/JHU 


Space Radiation Effects on Photodetectors 
for GAMES and GLAS 

Xiaoli Sun/JHU (rev. 2, Nov. 1993) 

SUMMARY 

GAMES and GLAS are expected to have a total radiation doses 
of 5 krad and 10 krad, respectively, with 100 mil thickness 
aluminum shielding, and 2 krad and 4 krad with 200 mil shielding. 
The peak dose rates are about 10* 3 rad/s with 100-200 mil shielding. 
The principle radiation particles are protons and electrons in the 
South Atlantic Anomaly (SAA) region and the cosmic rays near the 
two horns at the poles of the magnetosphere. The electron radiation 
can be attenuated to a negligible amount after 7 mm A1 shielding. 

The radiations of concern are the trapped protons in SAA and solar 
protons. 

The PMT in GAMES is expected to suffer little long term 
damage after 10 krad dose. The transient radiation induced noise 
counts are estimated to be 210 4 /s and 810 5 /s for a rad-hard and a 
non rad-hard PMT, respectively, which is not likely to cause 
significant degradation in the overall receiver performance since the 
dark counts due to stray light still dominates. However, radiation 
damages to the GaAs photocathods and microchannel plates still need 
to be studied. 

Si CCDs are inherently soft to radiation damages. The most 
significant damages are the reduction in the charge transfer 
efficiency and increases in the dark current due to surface ionization 
damages and bulk crystal lattice displacements. Some CCDs start to 
show noticeable performance degradation after less than 1 krad. 
Consequently, ordinary Si CCDs may not meet the requirement for 
GAMES, and therefore, rad-hard CCDs, such the CCDs in star trackers, 
should be used. 

Si APDs for GLAS may experience more than an order of 
magnitude increase in bulk leakage current, up to 720 pA for the 
1.064 pm wavelength enhanced APD and 60 pA for the 600 nm 
wavelength Geiger mode Si APD, after a total proton dose of 4 krad. 
The 1.064 pm altimeter receiver performance will be affected 
somewhat but the Geiger mode APD in the lidar receiver is expected 
to suffer series performance degradation due to the increase in dark 
counts, > 106 /s at the end of the five year mission. 
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I. Space Radiation Environment for GAMES and GLAS 
Spacecraft 

A. South Atlantic Anomaly (SAA), near 30° S. Latitude/345 0 E. 
Longitude (Messenger (1992), ch, 8), (Stassinopoulos (1988)). 
Mostly protons of energy 0.1-5 MeV and electrons of energy 
< 2 MeV. 

B. Cosmic Rays and solar flare protons near the two geomagnetic 
polar regions (the two horns at the poles of the magnetosphere). 
(Messenger, Stassinopoulos (1988), Barth). 

C. The total expected space radiation exposure for GAMES and GLAS 
missions are 5 krad/4yrs and 10 Krad/5yrs, respectively, with 
100 mils thickness shielding and 2 Krad/4yrs and 4 Krad/5yrs 
with 200 mils shielding. The total dose decreases very slowly as 
the shielding thickness increases to beyond 1/4 inch (500 mil). 
The peak dose rates are well below 10*3 rad/s with greater than 
100 mil shielding. The details are given by Barth and 
Stassinopoulos (1993a, 1993b). 

D. A review of radiation testing of semiconductor devices for space 
electronics has been given by R. L. Pease et al (1988). 

E. A general theory of radiation damages to microelectronics devices 
including opto-electronics devices have been given by 
Srour(1988) and Messenger (1992). 

II. Radiation Effects on PMT for GAMES 

A. Transient Effects 

1. Scintillation and fluorescence of the glass window, causing the 
dark current to increase proportionally to the radiation dose rate 
up to 10 1 ® rad/hr with recovery time of up to 30 ps ( Johnson, 
Viehmann, Levy). 
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2. Cerenkov emission, photons emitted by electrons traveling faster 
than light in glass (Birnbaum). Most of the Cerenkov electrons 
may be baffled out because they travel in specific direction and 
last only for picoseconds. The noise due to this effect appears as 
huge but short current spike at the PMT output and therefore can 
often be discriminated 

3. Secondary electron emission of the window material. 

4. Direct bombardments of electron radiation on the dynodes. 

5. Bremsstrahlung photons (photons released by the electrons 
decelerated in the shield). 

B. Long Term Damage 

1. No long term damage under low dosage (10 4 rads) (Johnson). 

2. Fluorescent and phosphorescence, lasting for days, weeks, and 
sometimes even months (Viehmann), causing the PMT dark 
current to increase by a few orders of magnitude (e.g. increased 
by about a factor of 10 after 10 5 rad (Johnson)). 

3. Minor changes in photocathod quantum efficiency (Johnson). 

4. Window brownish and higher transmission loss (Johnson). 

5. The photocathod itself suffers little damage because it is thin and 
has a low absorption coefficient to high energy particles. The 
damage to the dynode chain can also be neglected (Johnson). 

6. The radiation effects on PMTs are almost independent of the 
temperature (Johnson). 

C. Radiation Ha rdness. Shielding and Protection 

1. Minimizing the size and the thickness of the glass window. 

2. Proper choice of the window material (Johnson, Viehmann, and 
Birnbaum). 
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3. Minimizing the size of the photocathod to reduce the collection 
area of radiation induced photons. 

4. Electron focusing and baffling (Birnbaum). 

5. Anode current limiting to prevent bum-out by transit effect. 

6. Turning off the power at the peak of the radiation. 

7. Shielding with Aluminum, 6 mm thick would stop almost all the 
electrons and reduces the protons flux by more than one order of 
magnitude (Viehmann). More detailed data on shielding are given 
by Barth and Stassinopoulos. 

8. Very high energy particles and cosmic gamma rays cannot be 
shielded with reasonable amount of aluminum (Rasmussen). 
However, their occurrences are rare and only cause single event 
phenomena which is not a primary concern for GAMES and GLAS. 

D. Some Test Results of PMT for Galileo (Birnbaum et al (EMR)) 

1. Use of 60 Co Gamma Ray as test source which induces equal amount 

of florescence. 

2. The PMT is space qualified for the Galileo mission around Jupiter 
which has a harsher radiation environment than earth. 

3. The rad-hard PMT produced about 2 times 10? counts/s per rad/s, 
roughly 2.5% the noise count by a conventional 1" diameter PMT. 

E. Expected Radiation Damage to the GAMES PMT 

1. The radiation induced noise counts should be well below the 
estimated noise count due to stray light (lO^/s). 

2. The EMR PMT produce 10 3 * * 6 /s noise counts at about 0.4 rad(Si)/s. 

3. The peak dose rate likely to be encountered by GAMES is well 
below 10'3 rad(Si)/s with 100-200 mil shield according to Barth 
and Stassinopoulos (1993a, Fig. 25-27). The expected noise counts 
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due to radiation is therefore less than 20,000/s and 800,000/s for 
the EMR rad hard PMT and conventional PMTs, respectively. 

4. The GAMES receiver should be able to tolerate a long term dark 
count increase by more than a factor of 10 (a total dose of 10 3 * 5 rad 
according to Johnson). 

5. The GAMES PMT should suffer little radiation damage, even with 
non rad-hard PMT provided it is shielded with more than 

100 mils thick aluminum. 

F. New Issues to be Studied about Radiation Effect on GAMES PMT 

1. No test data was found on GaAs photocathods which are believed 
to be more sensitive to radiation damage. 

2. No test data was found about multichannel plates. 

3. There may be some test data about GaAs photocathod night vision 
instrument. 

III. Radiation Effects on CCDs for GAMES 

A, Transient Efforts 

1. Increasing noise counts due to ionization radiation, e.g. at 
4.3x1 0 J Vcm3 new hole-electron pairs per rad(Si) (Killiany, p. 

158). 

2. The statistics of the dark counts generally do not follow Poison 
distribution and have a much larger standard deviation (e.g. 30 
times for 1.25 MeV gamma rays) as compared with that of a 
Poison random variable (Killiany, p. 161). 

3. Well saturation, blooming into adjacent pixels, or even burn-out 
due to noise photocurrent spikes. Well saturations occur for dose 
rate less than 1.0 rad per pixel integration time (Killiany). The 

recovery time depends on the clock rate and may be as long as 

milliseconds (Killiany, p.159). 
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4. Single event phenomena, i.e. single pixel saturation or burn-out. 

B. Long Term Damages 

1. Ionization damages. Protons and electrons from the radiation 
create electron-hole pairs within the insulation silicon oxide 
under the gate, as in a MOSFET. Because of the different 
mobilities, electrons are swept away and holes left behind, 
causing a positive charge build up in the oxide. Radiation also 
causes increases in the interface state density (traps formed due 
to crystal lattice discontinuity between silicon channel and the 
gate oxide) (Messenger, p. 281) (Killiany, p. 156). 

a. Flat band (threshold) voltage shift 

1) 0.3 V/krad as reported by Roy (1989), 0.09 V/krad by 
Hopkinson (1991, 1992). About 0.01 V/krad for electron 
bombardment, as reported by Roy. 

2) There were no significant annealing effect on flat band shift 
(Hopkinson (1991, 1992)). 

3) The flat band voltage shift is much larger for positive gate than 
for negative gate voltages (n channel), because the charges are 
trapped near the gate oxide interface in the later case. 

4) The flat band voltage shift is larger at lower temperature due 
to the longer trap life time. 

5) The flat band shift was shown to be 3-4 times smaller if the 
devices were not powered (Hopkinson (1992)). 

6) Uniform flat band voltage shift up to a few voltages may be 
compensated by increasing the reset voltage. Non-uniformities 
in flat band voltage shift form traps and cause signal 
distortion (Killiany). 
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b. Charge transfer efficiency (CTE) decrease due to radiation induced 

interface state density increase. The CCD read noise will increases 

as CTE decreases. 

1) Charge transfer inefficiency (CTI=1-CTE) increased from 6xl0* 6 
to 3xl0* 5 after 340 rads of 6.5 MeV protons, from 10' 5 to 

8x1 O' 4 for high resistance devices after 2 krad of 10 MeV 
protons with little temperature dependence (Abbey). 

2) CTI dropped to twice the pre irradiation level after 16 hr 
140°C annealing (Abbey). 

3) CTI becomes too high for the CCD to operate after 10 5 rad for 
most of the applications (Killiany). 

4) Normal operation requires CTIxN«0.1 where N is the number 
of charge transfers (number of pixels in a column plus the 
number of shift register stages). 

5) Interface state trapping effects are much less for buried 
channel CCDs than for surface channel CCDs. 

c. Surface dark current increase due to interface state density 

increase (GAMES tracking CCD requires Idark«10 nA/cm 2 ). 

1) Dark current increase at 10 nA/cm 2 /Krad as reported by 
Hopkinson (1991, 1992). 

2) Pixel to pixel dark current non-uniformity was not significant 
(Hopkinson (1991)). 

3) The increase in dark current was about factor 2-4 less when 
the devices were not biased. 

4) The dark current increased by 12 nA/cm 2 after 100 krads of 
electron radiation (Roy). 

5) The dark current continued to increase after the radiation, the 
final dark currents became 2-3 times those immediately after 
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the radiation (Hopkinson (1991, 1992)). This antiannealing 
effect was not well understood and did not always occur (e.g. 
Abbey). 

6) The dark current decreased with temperature (Hopkinson 
(1992), Roy). 

d. A well designed CCD readout circuit suffered little radiation 
damage as compared with other parts of the CCD (Killiany). 

e. The power consumption of the CCD increased slightly with 
radiation dose (e.g. 0.04 to 0.2 mA/rad) (Hopkinson (1991)). 

f. It is widely assumed that equal doses of gamma rays and protons 

result in similar amounts of ionization damages, regardless of the 
proton energy (Raymond). Gamma rays are often used as test 
radiation sources. 

2. Crystal lattice displacement upon impacts by protons, forming 
defects which act as traps and charge recombination centers 
(Killiany, Hopkinson (1991, 1992)). However, the principle 
damages of CCDs due to protons and electrons are ionization 
damages. 

a. Dark current increases, 2.4 nA/cm 2 /MeV 3 months after 3 krads 
proton dose, with negligible pixel to pixel variation. Dark current 
due to bulk damage annealed after storage at room temperature 
(Hopkinson (1992,1992)). 

b. Charge transfer efficiency may decrease due to traps. The effects 
depend on the bulk trap emission time constant (0.53 ms as 
reported by Killiany) and the clock frequency. 

c. Degradation of minority carrier life time and concentration (no 
direct effect on imaging sensors but on the readout circuit. 

d. Protons only cause minor bulk displacement damages (which are 
the principle damages by neutrons). A dose of 1 rads protons of 
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1-100 MeV is equivalent to 1-2 times 10 7 n/cm 2 1 MeV neutron 
fluence (Raymond), 

e. The displacement damages due to protons starts to cause 
noticeable performance degradation at greater than 10 krads 
(about 10 9 * 11 n/cm 2 equivalent neutron fluence) (Killiany, p. 163). 

f. Displacement damage affects buried channel CCDs more than 

surface CCDs. 

g. Displacement damages due to Gamma rays are insignificant for 
dose less than 10^ rads(Si) (Killiany). 

C. CCD Radiation Hardness 

1. Aluminum shielding, 3 gm/cm 2 stops most of the electrons and 
reduces the protons flux by more than one order of magnitude 
(Stassinopoulos). 

2. Use of buried n-channel devices. 

3. Keep high uniformity across the entire imaging sensor array. 

4. Use of a threshold voltage tracker. 

5. Proper choice of the gate oxide. 

6. Turning off the power during the peak of radiation 
7 Current limiting to prevent pixel burn-out. 

8. Software correction or fault tolerance. Radiation induced noisy 
pixels appear scattered and instantaneous while the tracking 
beacon moves relatively slowly in a somewhat predictable way. 

9. In general. Si CCDs have initial observed degradation with 

potential failure after 10* rads, and severe degradation with high 

probability of total failure after 10 5 rads (Messenger, p. 642). 
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10. Modem rad-hard CCDs can be almost immune to ionization 
damage and its performance is primarily limited by the 
displacement damages (Dale et al). 

D. Expected Radiation Damage on the GAMES CCD Tracking Detector 

1. Transient effects: About 2,600 hole-electrons will be generated in 

a 20 by 20 by 15 pm 3 pixel per frame (0.1 seconds) at the peak of 
the radiation dose rate (10 -3 rad/s). 

2. Flat band shift: <0.1 V/Krad xlO Krad = 1.0V, which can be easily 

accommodated. 

3. Charge transfer inefficiency increases: CTI may drop from 10'5 to 
8xl0* 4 after 2 Krad, which is unacceptable. 

4. Dark current increases: lOnA/cm^/Krad x 4 = 40 nA/cm^, which is 
also unacceptable. 

5. The above results may be over pessimistic since most of the 
numbers quoted were drawn under the worst cases. 

Characteristics of radiation damages vary widely from devices to 
devices. A new test should be conducted for the GAMES CCD. 

6. The GAMES CCD and its shield have to be carefully designed or it 

will not meet the mission requirements. The primary damages 
are the increases in charge transfer efficiency and dark current 
due to ionization (0.5" shielding may have to be used to keep the 
radiation dose inside the shield well below 1 Krad). 

Ill Radiation Effects on Si APDs for GLAS 

A. Transient Effects 

1. Shot noise due to the radiation induced photocurrent, 

Irp= q go (G Vg + V s )(dy/dt) with q the electron charge, go the 
charge generation rate (4.04-1 0 13 /rad/cm 3 ), G the average APD 
gain, Vq the volume of the APD high field region, V s the volume 
of the APD which contribute to surface current, and dy/dt (rad/s) 
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the dose rate. For example, I rp =l 75 nA and I rp =156 nA under 
28 rad/s 60 Co gamma ray and 1.5 MeV electron radiation for a 
EG&G C30902S Si APD at G=100 (Swanson). 

2. Noise generation at the glass window, as in PMTs. 

B, Long Term Damages 

1. Increase in bulk leakage current due to displacement damages to 
the crystal lattice, at Albulk=a-4 ) - V, with a the leakage current 
constant, <}> the neutron radiation fluence, and V the volume of the 
active region of the APD (Kraner). 

a. Since 1 rad 1-100 MeV protons causes the same amount 
displacement damage as about 1.5'10 7 /cm 2 1 Mev neutron 
fluence (Raymond), Albulk=a- 1 .5- 1 0 7 -V (A/rad). 

b. Geiger mode Si APD (10V above the breakdown point), 
a=1.710-l6 A/cm (EG&G C30902S, V=610'6 cm3) (Buchinger), 
therefore, AIbulk=15 pA/krad. 

c. Analog mode Si APD (EG&G C30902S), Ibulk increased from 
0.01 pA to 3 pA after 200 Krad 1.5 Mev electron radiation 
(Swanson). However, the increase in Ibulk due to protons of the 
same dose should be several orders of magnitude greater 
(Van Lint, 1975). 

d. For the Si APD used in MOLA, d=800pm, L=140jim (Hammond), 

V=7.032 10’ 5 cm 3 , AIbulk=180 pA/krad for proton radiations. 

2. Increase in surface leakage current, primarily due to ionization 
damages at the APD surface between the insulation layer and the 
bulk (Swanson). 

a. The increases in I s should follow about the same rule as for Si 
MOSFET, in the order of 10 n A/cm 2 /krad. 

b. The surface leakage currents are known to exhibit 1/f noise 
(Swanson). 
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c. The noise due to surface leakage current is usually much smaller 
than other noise sources, e.g. bulk leakage current and 
preamplifier noise). 

3. Decreases in quantum efficiency due to surface leakage and carrier 
lifetime reduction due to radiations (Kalma, Wiczer, Barnes). 
However, this effect may be minimized by using a smaller load 
resistor for the photodiode. 

a. A reduction of quantum efficiency of 15% and 65% were observed 
after 10 4 5 * and 10 8 rad for Si PIN photodiode (Kalma). 

b. The lifetime reduction due to radiation displacement damage only 
affects low frequency response (Kalma). 

4. Damages to the preamplifier. GaAs preamplifiers use MESFET 
(JFET) and therefore are much more resistant to ionization 
radiation damages than Si MOSFET devices (Messenger). 

C. Photodiode Radiation Hardness 

1. APDs are much more sensitive to radiation damages than PIN 
photodiodes because Albulk is multiplied by the APD gain. Unlike 
PIN photodiodes, the APD dark noise is usually comparable with 
the preamplifier noise, and therefore, the increase in the APD 
dark current due to radiation are much likely to degradation in 
the receiver performance. 

2. The window material has to be chosen with care, as for PMTs. 

3. It is recommended that the shielding for the APD be increased to 
0.5" so that the net dose at the APD is below 1 krad for the GLAS 
mission. However, further increasing the shielding should have 
little effect in reducing the total dose (0.6 krad after 1.5"). 

4. Several rad-hard photodiode structures have been published 

(Wiczer, Barnes), but they may not be needed for the GLAS 


mission. 
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5. Si APDs were known to have passed radiation tests for free space 
laser communications, up to 10 5 rad. However those are all 
gamma ray tests done at EG&G (according to my conversation 
with Dr. McIntyre). Protons of the same dose are expected to 
cause much more long term damages to the APDs. 

D. Expected Radiation Damage to the Si APD in GLAS 

1. Total transient bulk dark current increase by about 0.056 pA 
under the expected 10- 3 rad/s peak dose rate, which is well 
below the average APD bulk dark current in analog mode. 

2. The transient increase in the noise counts in a Geiger mode APD 

for the lidar will be less than 0.056pA/q =355,000 counts/sec*^ 
which should not cause device saturation. ' 

/V4 . 

3. Permanent bulk dark current increase by about 0.72 nA and 

60 pA for the 1 pm enhanced Si APD and the Geiger mode APD 
for the lidar, respectively, after 4 krad total dose over the entire 
5 years GLAS mission (200 mil shielding). 

4. The increase in bulk leakage current will cause the APD dark noise 
spectral density to rise to a few times the preamplifier noise 

(2 pA/Hz 1 / 2 ) at the end of the mission, which will cause some but 
not catastrophic performance degradation. 

5. The radiation induced noise counts in Geiger mode APD in the lidar 
would greatly exceed the saturation level (10 7 counts/s) at the 
end of the mission regardless of the shielding thickness. 

6. Little loss in quantum efficiency at the end of the mission (4 krad 
total dose). 

7. The effect of increase in bulk leakage current on receiver 
performance is far more severe for APDs than for PIN 
photodiodes. The former is usually operated close to shot noise 
regime and the latter operated in preamplifier noise limited 
regime. 
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E. Some Radiation Test Data of the APPs for Laser Com and MOLA 

1. EG&G C30902E APD and preamp module by TRW was reported to 
have passed 10 5 rad(Si) gamma rays with no damage except a 
30% increase in the dark current (Conner). 

2. The specification of the APD for MOLA by McDonald Douglas said 

that the APD shall not be permanently damaged and shall 
perform within the specified limits after 10^ rad(Si) of radiation. 

3. The radiation dose should be assumed as the dose outside of the 
package. The APD itself inside the shield may have experienced 
less than 10 4 rad. 

F. New Issues to be Studied 


1. There has been no measurement data on long term damages due 
to proton radiations. Protons, which are the primary radiation in 
the GLAS orbit, are expected to cause significant increase in Ibulk 
due to displacement damages. Dr. McIntyre who is the leading 
expert in the field of APDs also expressed concerns about proton 
damages to Geiger mode APDs. 

2. It is recommended that we do a proton radiation test on Si APDs, 
especially Geiger mode APD photon counter. 
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1. Design of the GAMES PN Code Coarse Ranging Subsystem 

Figure 1 shows a block diagram of the GAMES coarse ranging 
subsystem design. The laser diode transmitter is modulated by both 
a sinusoidal signal as the subcarrier and a pseudo noise (PN) binary 
sequence as the subcarrier amplitude modulation. The former is for 
the fine ranging and the latter is for the coarse ranging. Figure 2 
shows the two possible combinations of these two modulation 
signals along with the hardware of how to realized them. Both of 
those modulation formats should result in the same coarse ranging 
performance. The depth of the PN code subcarrier modulation is no 
greater than 10%. 

Norman and Gardner [1] have shown that the optimal detection 
scheme for a PN code ranging system is the cross-correlation of the 
received photon counts with the PN code to find the maximum. One 
can also average the correlator output over a number of PN code 
period to improve the signal to noise ratio (SNR). To reduce the 
required computer time, one can equivalently form a histogram of 
the detected photon counts over a number of code periods and then 
cross correlate the histogram with the PN code, as shown in 
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Figure 1 . The cross correlation and the peak detection can be done in 
software. 

2a Coarse Ranging Performance 
2.1. The Ranging Resolution 

The resolution of the coarse ranging subsystem, Ad, is given by 

Ad=c/2x bin (i) 

where c is the speed of light and T bin is the duration of each bin of 
the histogrammer. The value of T bjn is usually equal to the bit time 
of the PN code. There may well be misalignment between the 
histogrammer and the received PN code. The largest timing offset is 
T b j n /2. This timing misalignment causes the peak value of the 
correlator output to decrease by as much as a factor of two for 
rectangular PN code pulse shapes. The amount of the noise of the 
correlator output is not affected. 

2.2. SNR at the Output of the Correlator 

The other important criterion of the ranging performance is the 
false alarm probability in locating the correlation peak. Because the 
low PN code modulation depth (sometimes called modulation index), 
the dominate noise source for the coarse ranging subsystem is the 
stray light and the unmodulated fine ranging signal. The statistics of 
the total noise in the correlator output can be assumed to follow the 
Gaussian distribution. The false alarm probability under this 
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assumption is solely determined by the SNR of the correlator output 
and the PN code length. 

The average detected photon counting rate can be written as 

^(t,T d )=i[(l-m PN )L p P s +P b ] +^[mp N LpP s .£_a,p(t-it-T d )] ( 2 ) 

where r\ is the photodetector quantum efficiency, hf is the photon 
energy, m PN is the modulation index, L p is the total signal round trip 
propagation loss, P s is the transmitted optical power averaged over 
a bit time, Pb is the background noise power, aj=0,1 is PN code, p(t) 
is the transmitted laser pulse shape, t is the PN code bit time, and 
Td is the time delay due to light propagation. The laser pulse shape 
is assumed to be a rectangle of unity amplitude and one bit time 
duration. 

The photodetector output is a Poisson random point process with 
the counting rate given by (2). The histogrammer output can be 
written as a random sequence, xo, xi,...,xn-i, with each element 
equal to the number of photons collected in that bin. It can be shown 
[2] that xo, x-|,...,xn-i are independent Poisson random variables with 
mean given by 

X k =Mx>.(kT,Td)T=M^mp N LpP s T a k+! + Mp^[(1-mp N )L p P s +P b ]T (3) 

where M is the number of PN code periods which the histogrammer 
has integrated and the propagation delay time, Td, is assumed to be 
equal to integer number of bit time, i.e. Td=lt. 

The signal output from the correlator can be written as 
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(4) 


where 


aj= 2a r 1 

is the bipolar form of the PN code. 


(5) 


Since the histogrammer output xo, x-i ,...,xn-i are independent Poisson 
random variables, the mean of the correlator output can be derived 
as 


t *iaLj 


i = 0 


=M^m PN L p P s TZ a k+1 a| +j + M^Kl-m^LpP.+PjT Z a| +j 


(6) 


When a maximum PN binary sequence of code length N is used, it has 
the following property 



(7a) 


and 



j =0,1,2,... 
else 


(7b) 


The mean of the correlator output becomes 

9 _/ M H7f C 1 -m PN > L p p s + p bjT + Mjlm^LpP.fN+l, j = | 

yj ^M^(1- mpN )L p P s+ P b]T , jw (8) 
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The useful signal is defined as peak value of the correlator output, 
yj»i, minus the constant background, i.e., 

y s =lM(N+1)im PN L p P s x (9) 

The variance of the correlator output can be calculated as follows. 

Since the histogrammer output xo, xi xn-i are independent Poisson 

random variables, 

VARlYj !=X (a| + j ) 2 -VAR{Xi}=1:’ VARIxiIkU X, (10) 

1 = 0 1 = 0 i = 0 

Substituting (3) into (10) 

VARjyj }=Mim PN L p P s T N f a ltl + MNif(1-m PN )L p P s+ P b ]x (11) 

Since X , 

VARIyj l=c y 2 =lM(N + 1)im PN L p P s T + MN^KI-m^LpP.+PbjT 

( 12 ) 


Notice that the variance is independent of the bin number. 

A signal to noise ratio can be defined as 

SNR=|j (13) 

Substituting (9) and (12) into (13) 
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SNR =< 


M(N+1) 

2 


mp^LP s 


TIT 

mp Nhf*-p^s + j^$0-"WL p P s +P b ] 


2N Wr 


(14) 


It is noticed that (14) is different from that given in [4] which we 
believe is incorrect. 


If we define the average number of detected signal and background 
noise photons per bit time as 


"s=^ L pP s 


(15) 


and 


n -Hip 
b hr b 


(16) 


the correlator output SNR becomes 

SNR =n /- M ( N 2 +1 > . 


rnp^n. 


\/ m PN n s + ~ m PN) n s +n b] 


( 17 ) 


2.3. SNR for Non Photon Counting Receivers 

The photodetectors, either a PMT or an APD, can operate in analog 
mode. The photodetector output in this case is simply the 
photocurrent integrated over a bit time. The receiver does not 
contain a discriminator but an integrator or a filter which perform 
the integration. Besides the shot noise associated with photon 


6 


February 1 7, 1 994 


X. SurVJHU 


detection, there is also the so called excess noise which results 
from the randomness of the photodetector gain. When an APD is used, 
the preamplifier noise also has to be considered. 

The mean of the histog rammer output when an APD is used can be 
written similarly to Eq. (3), as 

s|H^GX(kT,T d )T + ^ = 

= MG^m PN L p P s T a k+ | + MG^(1-m PN )L p P s+ P b ]x + MG^ + 

(18) 

where G is the average photodetector multiplication gain, t| a is the 
quantum efficiency, and Ids and Idb are the photodetector surface and 
bulk leakage current, and q is the electron charge. Notice that the 
quantum efficiency of an APD, T| a , is normally a few times photon 
counting probability in Geiger mode. The variance of the 
histogrammer output can be written as [5] 

Var{x k }=MG 2 F|^mp N LpP s T a k+ | + -^j[(1-nip N )LpP s +P b ]T + 

+ (19) 

where N am p is the one sided spectral current density of the 
preamplifier noise in A 2 /Hz and F is the APD excess noise factor 
given by [6] 

F ~ k e ffG+(2-1/G)(1-k eff ) (20) 
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with keff the effective APD hole to electron ionization coefficient 
ratio. 


Substituting (18) and (19) into (6) and (10) and following the 
procedure given in the previous subsection, the correlator output 
SNR can be written as 

SNR' =N y / M ffl+ 1 ) 


m 


PN* 


— L P 

*-p r s 


hf 


aT [m 


hf 


PN 


L p P s 


1 -m PN )L p P s+ P b ] + Flf Jj 




N T 
2q 2 ' 
(21) 


|2 


or 

SNR' =sj M (N+ 1 ) . 

m PN n s 

\/ ^pnHs ^[(1 -m™K+n b+ ^]] + + 

( 22 ) 

When a PMT is used, Eq. (22) can still be used to calculate the 
correlator output SNR provided an appropriate excess noise factor F 
(F-1.1-1 .5) is used. The gain of a PMT is usually so large (>10 5 ) that 
the preamplifier noise can be ignored. 


For base band PN code ranging system, there is no subcarrier. The 
laser diode is directly intensity modulated by the PN code. One can 
still use Eq. (17) and (22) to calculate the correlator output SNR 
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provided the modulation index, mpN, is replaced by the laser OFF to 
ON extinction ratio. 


2.4. Probability of Incorrect Measurement vs. SNR 


The probability of the incorrect range measurement can be written 
as 


P E = Prob(yi < Y;, 0< i < N-1 , i * I) (23) 

As an approximation, we can assume that the yi's follow the 
Gaussian distribution with the mean and variance given by (6) and 
(12). The probability of the incorrect range measurement can be 
written as 


P P =1- 


‘/ZKGy 


-exp[ ^ ^ 


r 


2o^ 




•/2jca, 


2 

-exp[^%^-]dy 2 


2a£ 


N— 1 


dyi 


(24) 


or 


-'-ir 


e-y 


1 -ierfc (y+^jr) 


N-1 


where erfc(x) is the complimentary error function defined as 

erfc (x) = 1 - erf(x)=1-J JLe _vZ dv=£ ^-e _v2 dv 

For large value of SNR, 

P (E - e-y 2 erfc (y+^)dy . 


(25) 


(26) 


(27) 
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3. Numerical Results 

The probability of ranging error was evaluated numerically using the 
software called Mathcad and the results are given in the attached 
computer printouts. 

The first attachement shows the data and the plot of Pe vs. SNR for 
code length N=29-1 =51 1 . It shows that a SNR of 18 dB is required to 
achieve Pe < 10* 6 for N=51 1. This result can be scaled for other code 
lengths because Pe is proportional to the code length for any given 
SNR, as shown in (27). 

The second attachement shows the calculated performance of the 
GAMES PN code coarse ranging sybsystem. A photon counting receiver 
is assumed. All the system parameter values used in the calculation 
are given in the printout. The total histogrammer integration time is 
taken to be 0.1 second. If we assume the target is moving at 1 .0 m/s, 
the timing offset between of the histogrammer and the received PN 
code due to the Doppler shift will change by 2/3 ns (3.3% of the bit 
time) during the 0.1 second integration time interval. Furthermore, 
the SNR is reduced by one half in the calculation to account for the 
possible misalignment between origins of the histogrammer and the 
received PN code. The calculation results show that a detected 
signal photon rate of 1.9-106 is required to achieve Pe < 10-®. The 

dominate noise sources are the shot noises from the signal and 
background light. 
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The third attachment is the same as the second but the laser 
modulation index was changed from 10% to 5%. A detected signal 
photon rate of 5.8-1 0 6 is required to achieve Pe < 10* 6 in this case. 

The fourth attachement shows Pe vs. detected signal photons per bit 
time for a base band PN code ranging lidar system which use an APD 
photodetector in analog mode. The system parameter values used are 
typical for practical devices available today. The dominate noise 
source in this case is the preamplifier noise. The calculation results 
show that 2 photons/bit time is required to achieve Pe < 10' 6 for 
code length N=2 7 -1=127. 

For the same histogrammer integration time, the SNR at the 
correlator output is virtually independent of the code length 
according to Eq. (22). On the other hand, the ranging error probability 
is proportional to the code length for any given SNR. The code length 
should be determined by the maximum range the target or the 
minimum range unambiguity required by the system. 
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2GHz + simulated Doppler shift 



Figure 1 . Design of the breadboard of the GAMES coarse ranging subsystm. 
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Figure 2. PN code subcarrier modulation for the GAMES coarse ranging subsystem. 


Probability of Error vs. SNR 
of a PN Code Laser Ranging System 

N := 51 1 The code length 
SNR := 5, 5.5.. 15 


PE( SNR, N) := 


N- 1 
2 (n 


20 


exp(->' 


•20 


,’ 2 )- |l - erf^y + 


SNR 


r* 


w * 


TOL := 10 


15 


The integration tolerance 


SNR 20- log( SNR ) PE(SNR,N) 


0.1038 

0.0257 

5.6331- 10" 3 

1.0972- 10" 3 

1.8949-10“ 4 

2.9-10" 5 

3.9314-10“ 6 

4.719- 10~ 7 

5.0137- 10“ 8 

4.7137- 10“ 9 

3.9205-10“ 10 

2.8841-10“ 11 

1.8762-10“ 12 

1.3241- 10“ 13 

5.6189-10“ 15 

0 

0 

0 

0 

0 

0 


5 

5.5 

6 

6.5 

7 

7.5 

8 

8.5 

9 

9.5 

10 

10.5 

11 

11.5 

12 

12.5 

13 

13.5 

14 

14.5 


13.979 

14.807 

15.563 

16.258 

16.902 

17.501 

18.062 

18.588 

19.085 

19.554 

20 

20.424 

20.828 

21.214 

21.584 

21.938 

22.279 

22.607 

22.923 

23.227 

23.522 


WRITEPRN( sradai ) : = 20- log( SNR ) 
WRITEPRN( pedal ) :=PE(SNR,N) 







GAMES PN Code Coarse Ranging Subsystem Performance 


PN CODE PARAMETERS 
Code Length: M 0ftl7 

T := 0.1 
M = 2.442* 10 3 

DETECTED SIGNAL AND NOISE PHOTONS 

Detected signal photons per bit time: 

n s := 0.02, 0.0225.. 0.05 


Code bit rate (bps): BR := 50- 10 6 

1 — s 

Bit time (second): z := — t = 2*10 

BR 

Total Integration time of the histogrammer (second): 
Number of code length averaged: 

M 1= floor 



Laser modulation index: m pN := 0.10 

Detected background noise photons per bit time: 

n b := 10 6 -t (assuming 10 A 6/sec) 

n b = 0.02 

PMT PARAMETERS 

Dari< count per bit time: n d :=6.3-10 4 -t 

n d = 0.00126 



CALCULATION OF THE NOISE VARIANCES 


Shot noise due to the signal: 


N s (n s ) := 


mpN ’ n s + 7?7 (1 ~ m PN ) ' n s 


2-N 


N .( 0.1 ) = 0.18991 


Shot noise due to background light: 

M - 2 ' N 

b N+ 1 " b 


Nu = 0.03998 


Shot noise due to dark counts: 


N j := 


2-N 


d -]7Tr n d 


N d = 0.00252 


CORRELATOR OUTPUT SNR: 


V 


SNR( n .) := 


M-(N+ 1) 

2 m PN' n s 


^N s (n s ) + N b +N ( 


SNR(O.l) = 32.80153 


Considering the loss due to misalignment between 
the histogrammer and the received PN code: 


SNR j(n s ) :=--SNR(n s ) 


CALCULATION OF RANGING ERROR PROBABILITY: 


PE(n ,) •'= 


N- 1 

2 fn 


20 


-20 


exp(-y 2 ) 




dy 


TOL := 10 


15 


The integration error tolerance 


NUMERICAL RESULTS 


T 10 6 



0.02 


0.0225 


0.025 


0.0275 


0.03 


0.0325 


0.035 






[0.0425 




0.0475 


0.05 



SNR j(n s ) 


5.57409 


6.09369 


6.58971 


7.06469 


7.52076 


7.95975 


8.38321 


8.79251 


9.18881 






201og(SNR ,(n s )) PE(n s ) 


14.92349 


15.69761 


16.37733 


16.98186 


17.52523 


798 



18.46821 


18.88226 


19.26519 


19.62111 


19.95337 


0.08284 


0.01678 






0 0.25 0.5 0.75 1 1.25 1.5 1.75 2 2.252.5 


x 10 6 


WRITEPRN( nsdai ) :=■ 


WRITEPRN( pcvsns ) := PE( n ) 





















































GAMES PN Code Coarse Ranging Subsystem Performance 


PN CODE PARAMETERS 


Code Length: 


N := 2047 


Code bit rate (bps): BR := 50- 10 6 

1 — C 

Bit time (second): x := — x = 2*10 

BR 

Total Integration time of the histogrammer (second): 
Number of code length averaged: 

M floor 



T := o.i 
M = 2.442* 10 3 


DETECTED SIGNAL AND NOISE PHOTONS 

Detected signal photons per bit time: 

n s := 0.05, 0.06.. 0.2 

Laser modulation index: m pN := 0.050 

Detected background noise photons per bit time: 

n b := l0 6 -x (assuming lO A 6/sec) 

n b = 0.02 


PMT PARAMETERS 


n d :=6.310 4 x 
n d = 0.00126 


Dark count per bit time: 



CALCULATION OF THE NOISE VARIANCES 


Shot noise due to the signal: 


N s (n s ) := 


2-N 


m 


p N’ n s + 1 ” m PN^‘ n s 


N-(O.l) = 0.19491 


Shot noise due to background light: 

XI - 2 ~ N 

b '~N+ 1 " b 


N b = 0.03998 


Shot noise due to dark counts: 


XI - 2 ' N 

d N+ 1 n d 


N d = 0.00252 


CORRELATOR OUTPUT SNR: 


SNR(n-) := ' 


M-(N+ 1) 


•m pN -n s 


I N s (n s ) + N b +N ( 


SNR(O.l) = 16.22731 


Considering the loss due to misalignment between 
the histogrammer and the received PN code: 


SNR j(n s ) :=--SNR(n s ) 


CALCULATION OF RANGING ERROR PROBABILITY: 


PE(n s ) := 


N- 1 

2- / "n 


20 


1-20 


exp(“y 2 )- [l - erfly + 


SNR 

~7i 


=1 


dy 


TOL := 10 


15 


The integration error tolerance 



NUMERICAL RESULTS 


Hl.J. 


n s T 10 6 SNR j( n s 


0.05 


2.5 


5.28375 

0.06 


3 


5.94033 

0.07 


3.5 


6.54205 

0.08 


4 


7.09994 

0.09 


4.5 


7.62186 

0.1 


5 


8.11366 

0.11 


5.5 


8.57977 

0.12 


6 


9.0237 

0.13 


6.5 


9.44821 

0.14 


7 


9.85556 

0.15 


7.5 


10.24761 

0.16 


8 


10.62591 

0.17 


8.5 


10.99178 

0.18 


9 


11.34633 

0.19 


9.5 


11.69054 

0.2 


10 


12.02523 


201og( SN R j(n s )) PE(n s ) 


14.45884 

15.47621 

16.31428 

17.02509 

17.64122 

18.18433 

18.66952 

19.10769 

19.50699 

19.87362 

20.21245 

20.52732 

20.82136 

21.09711 

21.35669 

21.60187 


0.19114 

0.02725 

0.00381 

5.27459-10" 4 

7.22913- 10" 5 

9.84565-10" 6 

1.33472- 10‘ 6 

1.80314-10" 7 

2.42957-10" 8 

3.26706-10" 9 

4.38643-10" 10 

5.88226-10" 11 

7.88051-10" 12 

1.05519-10“ 12 

1.62124-10" 13 


1.91544- 10" 14 



WRITEPRN( nsdat ) := S — WRITEPRN( pcvsns ) := PE( n ) 

f 10 6 S 








Probability of Ranging Error vs Detected Signal Photons per Bit 


PN CODE PARAMETERS 
Code Length: n:=127 

Code bit rate (bps) : BR : = 1 . 5- 1 0 6 

Bit time (second): x :=— x = 6.66667*10 7 

BR 

Total Integration time of the histogrammer (second): T 1= 0.1 

Number of code length averaged: 

M := floor^-Q M = 1.181* 10 3 


DETECTED SIGNAL AND NOISE PHOTONS 

Detected signal photons per bit time: 
n s ;= 1,1.25.. 3 

Laser modulation index: m pN 1= 0.85 

Delected background noise photons per bit time: 
n b := 10 

APD PARAMETERS 

APD Gain: G :=200 

APD ionization coefficient ratio: 

APD excess noise factor: 

APD bulk leakage current (A): 

APD surface leakage current (A): 

PREAMPLIFIER NOISE 

Preamplifier spectral noise current density (A A 2/Hz): 

(Analog Model 311 precision low noise current amplifiers) 

N iamp := ( I-O'IO - ' 2 ) 2 


k := 0.02 

F^k^G+^-lj-d-keff) 
F = 5.9551 

l(fc := i -0- 10~ 1 2 

Ids := 151°“ 9 



CALCULATION OF THE NOISE VARIANCES 


Shot noise due to the signal: 


N s (n s ) :=F- 


2-N 1 

m PN' n s + 1 "" m PN^ n si 


N.( 10) = 68.34408 


Shot noise due to background light: 

vj 2-N 

b ‘ F N+ 1 " b 


N b = 118.17152 


APD dark current noise: 
Electron charge (C): 


q := 1.6-10 


ri9 


N d :=JiL.( F .i*2 + !*Z] 

N+1 l <1 ,o 2 J 


N d = 52.33872 


Preamplifier noise: 


XT 2-N N iamp 
N a -TT^- 


N+ 1 2 . G 2 q 2 


N a = 645.9554 


CORRELATOR OUTPUT SNR: 


SNR(n,) := 


_ 1 


M-(N + 1) 

2 m PN' n s 


J N s (n s ) + N b +N d +N ; 


condering the loss due to misalignment of the histogramer 

1 
2 


SNR !(n s ) :=i-SNR(n s ) 


CALCULATION OF RANGING ERROR PROBABILITY: 


PE(n s ) := 


N- 1 

*r* 


20 


“20 


exp(-y 2 ) 


H'lr 1 )} 


dy 


TOL := io~ 15 


The integration error tolerance 


NUMERICAL RESULTS 


SNR j( n $ ) 20-log(SNR ,(n )) 


4.07216 I | 12.1965 


14.12569 




17.0303 


18.18119 


19.19531 






21.66741 


PE(n $ ) 


0.25097 


0.02039 





WRJTEPRN( pcvsns ) PE( n $ ) 
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Power Spectrum of PN Ranging Codes 

Xiaoli Sun/JHU, (Rev. 7 ) April 7 994 

A PN code sequence used in a ranging instrument can be written as 

x(t)= Z a k p(t-kAt) 

k = -oo 

where ak = ± 1 is the value of the k th bit of the PN code sequence, p(t) is 
the pulse shape, and At is the bit interval. Here we assumed that the PN 
code is M bits long and repetitive. 

The power spectrum of PN code sequence is given by [1 . p. 35] 

S(cd)=J R(i)e‘-i WT dT 

where R(t) is the autocorrelation function of the PN code defined as 

R(i)=lf T/2 x(t)x(t+i)dt 

1 J-T/2 

1 . Power Spectrum of Maximum PN Ranging Codes 

When a maximum PN code [2] is used, the autocorrelation function can be 
written as 

R " M =T n ?, <^- nT > - R 

where 

p (t)=^ J p(t)p(t+i)dt 

Substituting (4) and (5) into (2) 

s ^=^ lp (<<L e ' jn “ T -f «(<•>) 
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where P(co) is the Fourier transform of the pulse shape, i.e., 

|P(w)| 2 =J p(i)e"- il0T dT=| I* p(t)e _J <ot dt 


Since [3] 


n e-J noT = f ~ f n 5(t-nT)]e"j ot dt = y L n E_ 5(a)- nip) 


(7) 


( 8 ) 


the power spectrum of a maximum PN code, Eq. (6), becomes 

S(a,) = ? ^t2 ) |P(m)|; n JL n ¥ } " Tf ' 5(0l) < 9 > 

Therefore, the power spectrum of a maximum PN code consists of a series 
of discrete components at multiples of the repetition frequency of the PN 
sequence. The envelop of the spectrum is equal to the power spectrum of 
the pulse shape, IP(co)l 2 . 

2. Power Spectrum of JPL Ranging Code 

The JPL ranging code [4] is formed by modulo 2 sum of several shorter 
binary code sequences called components. It has the advantage of fast 
target acquisition. 

Let us first consider a JPL code, Z={zo, z-],...,zm}, which has two component 

code sequences, X={xo, xi xmi} and Y={yo, yi yM2h i.e., Z=X©Y. Here Xj, 

yi, and Zj are binary variables and equal to either 1 or zero. The length of 
the combined code is M=M1- M2 where Ml and M2 have to be relative prime 
numbers. 

The autocorrelation function of the binary sequence Z is given as 
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Rz(n)= R i?o Zi ® 2|+n= 


MT^I 

Ml xM2 ii? 0 


1 


M2^1 

j2 I o (Xj 1 ey i2 )©(x il+n ©y i2+n ) 


MT^l 

Ml xM2 ii? 0 


1 


M2^1 

j2 Z o (x il ©x jl+n )©(y i2 ©y j2+n ) 


( 10 ) 


Eq. (10) may also be interpreted as the difference in the probabilities of 
"1" and "0" for the binary variable Zj~z, +n , or, equivalently, 
(Xii0X il+n )e(yi2@Xi2 +n)» 2S, 


R z (n)=p(0)-p(1 ) (ii) 

Similar to the example given in [1 , p. 75, Figure 5.8], a Karnaugh chart can 
be used to help to calculate those probabilities. 


Yi2, yi2-n 


( X 1 i©Xji +n )©(yi 2 ©Xj2+n) 

00 

01 

1 1 

10 


00 

0 

1 

0 

1 

Xii, xjT+n 

01 

1 

0 

1 

0 


1 1 

0 

1 

0 

1 


10 

1 

0 

1 

0 


The binary sequence X and Y can be considered as statistically 
independent, i.e., p(Xn,X nw ,yi 2 ,y a+n )=p(x n ,Xi 1+n )-p(y a ,y i2+n ). For two 
relatively long component codes, the probability for (xn, x,i +n ) or (yi 2 , 
yi2-n) to have each of the four combinations are approximately equal and 
so is the probabilities of each outcome of the binary variable 
(xij©x il+ n)e(yi2©Xi2+n). 
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When n is a multiple of the code length of either, but not both, component, 
two of the columns or rows of the above karnaugh chart are deleted. 
Nevertheless, the probabilities of each outcome for 
(xi j©xn +n)©(yi 2@ x i2+n) are still approximately equal. When n is a 
multiple of both component code lengths, the value of 
(xTjSxn_ n )S(yi2SXi2-rn) ' s always unity. 


Therefore, the autocorrelation function of the JPL ranging code sequences 
can be approximated as 


R z (nH n 

' 0, otherwise 


( 12 ) 


The error caused by the approximation decreases as the component code 
lengths increase. 

The above result also applies to JPL code of arbitrary number of 
components. One can consider any one component and the combined 
sequence of the rest components as the two components. The two code 
lengths are still relative prime and each codes still consists of 
approximately equal numbers of 0's and 1 's. 

Since the autocorrelation function of a JPL ranging code is the same as 
that of a regular maximum PN code, their power spectrums must be the 
same. 


It is the crosscorrelation functions in which JPL ranging codes differ 
from regular maximum PN ranging code. 
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